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Abstract. Certain Bernoulli convolution measures /i are known 
to be spectral. Recently, much work has concentrated on determin- 
ing conditions under which orthonormal Fourier bases (i.e. spectral 
bases) exist. For a fixed measure known to be spectral, the ONB 
need not be unique; indeed, there are often families of such spectral 
bases. 

Let A = J- for a natural number n and consider the Bernoulli 

In 

measure with scale factor A. It is known that L 2 ([i\) has a Fourier 
basis. We first show that there are Cuntz operators acting on this 
Hilbert space which create an orthogonal decomposition, thereby 
offering powerful algorithms for computations for Fourier expan- 
sions. 

When L 2 (/j,\) has more than one Fourier basis, there are natural 
unitary operators U , indexed by a subset of odd scaling factors p; 
each U is defined by mapping one ONB to another. We show 
that the unitary operator U can also be orthogonally decomposed 
according to the Cuntz relations. Moreover, this operator-fractal 
U exhibits its own self-similarity. 
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1. Introduction 



In the theory of fractals and fractal measures, patterns of symme- 
try and self-similarity abound. In this paper, we study a new class of 
unitary operators on fractal L 2 spaces. This class of operators have sur- 
prising self-similar properties which, to our knowledge, have not been 
observed before. In one special case, we discover even more intricate 
patterns in the matrix of the operator. 

Fractal sets which are invariant under a collection of contractive 
maps (an iterated function system) exhibit scaling self-similarity; one 
sees the same shape when zooming in to look more closely at one part of 
the set. We might call this scaling in the small to find self-similarity. It 
is also common to find that the discrete sets which index Fourier bases 
for fractal L 2 spaces have an expansive self- similarity themselves. We 
say that these sets (called spectra) have a self-similarity in the large. 
When a measure \x has such a discrete index set V for a Fourier basis, 
we call (/i, T) a spectral pair. 

For illustration, consider first the simplest case — the Bernoulli con- 
volution formed by recursive scaling by ~ with two affine maps on the 
real line. The resulting measure /xi is an infinite Bernoulli convolution, 

4 

also called a Cantor measure, or a Hutchinson measure. The Hilbert 
space L 2 (/ii) has a Fourier basis which has a self-similarity under scal- 
ing in the large by 4. It is somewhat surprising to find that that ui 

4 

also gives rise to a symmetry based on scaling by 5. It turns out that 
scaling by 5 transforms the Fourier spectrum T into another ONB 5r. 
As a result, we have a natural unitary operator U acting in L 2 (/ii) 
which maps one ONB to the other. Its spectral properties turn out to 
reveal a surprising level of symmetry and self-similarity which lead us 
to the nomenclature operator-fractal. 

After briefly reviewing Bernoulli convolutions and some of their prop- 
erties in Section [2J we demonstrate the intimate connection between 
the spectral theory of affine fractals and representations of the Cuntz 
algebras from the theory of C*-algebras in Section [31 In particular, 
the scaling operations on the spectrum V influence the structure of the 
Hilbert space L 2 (n). We find that the scaling operators are isometries 
which satisfy Cuntz relations (Definition I3.ip ; the Cuntz isometries in- 
duce the orthogonal decomposition of L 2 (jj) shown in Proposition 13.51 
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For certain odd integers p, the scaled set pT also is a spectrum for a 
Fourier basis on L 2 (/i) (see |LW02l IDHS091 IDJ091 IJKS11] for details). 
We denote by U the unitary operator which enacts the scaling by p on 
the ONB r. In Section [4j we examine the interplay between U and 
the isometries that form the Cuntz algebra. In Section 14. 2\ we prove 
commutation relations for U and one of the Cuntz operators Sq. More 
applications of Cuntz algebras to the setting of fractal measures can 
be found in [BJ99] . 

In our main theorem, Theorem 14. S\ we show that when (^-,pT) is a 
spectral pair, then the associated operator U is an orthogonal sum of 
a rank-one projection and infinitely many copies of the same operator; 
these copies arise from the Cuntz decomposition. In the special (4, 5) 
case, the operator which appears infinitely many times in U is just a 
unitary multiplication operator on U itself (Theorem 14.10)) . Finally, 
we study the structure of the matrix of U acting on the first Cuntz 
subspace in the (4,5) case (Theorem 14.151) . 



2. Background and Notation 

In this section, we outline the construction of the Bernoulli convolu- 
tion measures and provide some of the necessary background results. 

2.1. Affine iterated function systems. The Bernoulli convolution 
measures are a special class of fractal measures, those generated on 
the real line by two affine maps with scale factor A G (0, 1). There 
is a rich literature on these measures; see |Erd39[ IPS96) IPSSOOl ISid03 



ILW021 IJKS071 IJKS081 IDJ091 IDHS091 IDHSWlll [JKSTT] for just a small 
sampling. The problem we consider here fits into a wider framework of 
spectral pairs, self-similarity, and spectral duality. For the interested 
reader, other related work in the literature includes Fug74, ITao04| 



ICN981 ICY061 IPW0H IWan02l iLiOTj . The interplay between the scale 
factor A and a spectral scaling between orthonormal bases has been 
considered in |JR95j . 

A Bernoulli convolution measure p,\ on the real line can be realized 
as an invariant measure for an iterated function system (IFS). Let 
{r + ,r_} be the IFS 

(2.1) t+(x) = X(x + 1), r_(z) = A(z-l), 

where the scaling factor A is in (0, 1). 

By the well-established algorithm from |Hut81j . the IFS {t + ,t_} 
generates a compact attractor set, X\, that satisfies 

X A = r + (X A )Ur_(X A ). 
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The IFS also generates a measure fi\ (the Bernoulli convolution), sup- 
ported on X\, which satisfies a similar invariance equation: 



(2-2) ^ x= 2 v A ° T+ / + 2 V A ° T_ 

We might ask under which conditions the measure n\ has a Fourier 
basis (i.e. an orthonormal basis of complex exponential functions) for 
the Hilbert space L 2 ([i\). When such an orthonormal Fourier basis 
exists, we say that fi\ is a spectral measure. 

This line of questioning has its origins in |JP98j , in which Jorgensen 
and Pedersen demonstrated Fourier bases for L 2 (/ia) when A = w~ 
for each n G N. They also showed that when A = ^+1' there is 
no orthonormal basis (ONB) consisting of exponential functions, and 
in fact, every orthogonal collection of exponentials is finite when the 
denominator of A is odd. The Fourier basis found in |JP98l Cor. 5.9] is 
indexed by a discrete set r C R given by 



(2.3) r = r (^J = \l^ a ^ 2n Y : e |0,-^,m finite 



For example, when n = 2, i.e. A = |, the spectrum is 

if" 1 1 
r (4) = | E a ^ : Oi e {0, 1} | = {0, 1, 4, 5, 16, 17, 20, . . .}. 

Note that the elements of V are integers when n is even and are all in 
|Z when n is odd. 

Remark 2.1. In the later sections, we will be considering other ONBs 
for the same Hilbert space L 2 (fi). To distinguish the original Jorgensen- 
Pedersen basis, we will call T the canonical ONB for L 2 (n) and label 
the others alternate ONBs. 

We observe that T has a self-similarity by scaling in the large. For 
instance, T is invariant under scaling by the value 2n, the reciprocal of 
A. But there is even a stronger scaling invariance: 

2nT U (2nT + ™ 



Recent results have shown that if A > |, there is no Fourier basis 
[DHJ09t Thm. 3.4]. It has also been shown that if A = j- where q is 
odd, the set of exponentials shown in |JP98j to be an ONB for 
are also orthogonal for n_a_ jJKS08l Thm. 11.9], |HL08l Thm. 1.2]. 
There are still a variety of open questions regarding the existence and 
classification of Fourier bases for Bernoulli convolution measures. 
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We note that, since we wish to discuss ONBs for Bernoulli measures, 
we will be restricting here to the case where the scale factor A = j- for 
n6N. To keep the notation simple, we will write \x for fij_ and X for 
Xj_ . We will use the notation e 7 to represent the exponential function 
e 27"7(-)_ if a collection of exponential functions {e 7 : 7 G T} indexed 
by the set T forms an orthonormal basis for some L 2 (fi\), we call T a 
spectrum for L 2 (fi\). Also, for ease of notation, we will write E(F) for 
the set of exponentials {e 7 : 7 G T}. 

2.2. The Fourier transform /i and its zero set. Given the invari- 
ance equation (12.21) . there is a standard convenient expression for the 
integral of an exponential function e t . We denote the resulting function 
in t by Ji since this produces a Fourier transform of fi: 



j2(t) = [ e 2nixt dfi(x 



x x 

gMACx+l)* + e 2ni\{x-l)t \ d ^ 



A 



1 

2 

= cos(27rAt)/i(At) 

= cos(27rAi)cos(27rA 2 t)/i(A 2 t) 



Continuing the iteration, we find an infinite product formula for Ji: 

00 

(2.4) /I(t) = JJcos(27rA fe t). 

k=l 

Given exponential functions e 7 and ey, we note that 
(ey,e 7 ) L2(M) = / e 2 ™^'> d/x(x) 

(2.5) = £( 7 -Y). 

We are considering orthogonal collections of exponential functions, so 
we consider the zeroes of the function /I By Equation (12.41) . /2 is zero 
if and only if one of the factors in the infinite product is zero. (This is 
not difficult to show, but see |JKS08j Lem. 11.2] for explicit details.) 
The cosine function is zero at the odd multiples of |, which yields the 
set of zeroes for /2, denoted Z: 



(2.6) Eft 



(2n) fc (2m+ 1) 



m G Z, k > 1 
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Given a discrete set T, then, the collection of exponential functions 
E(T) is an orthogonal collection if (ey , e^) L 2^ = 5 7i y, which occurs if 
and only if 

7 — 7' G Z for all 7, 7' G T with 7 7^ 7'. 

Note that the set from Equation ( 12.31) does indeed satisfy this condition. 
We computed earlier that 

£(t) = cos(27rAt)/i(At). 

Making use of our current situation for which A = and a change of 
variables, we also have the useful equation: 

(2.7) fi(2nt) = cos(27rf) fi(t) 

for all tel. The utility of this formula comes when t is an integer, as 
then we can actually cancel the factor 2n, i.e. 

jl(2nt) = pi(t). 

3. Cuntz Operators 

In this section, we define operators on L 2 (fi) which satisfy Cuntz 
relations. The resulting Proposition 13.51 will provide an orthogonal 
decomposition of the Hilbert space L 2 (fi), which allows us to study 
Fourier bases for L 2 ([i) in Section HI In many cases, there may be more 
than one spectrum for L 2 (fi), and if so, we can define quite naturally the 
unitary operator U mapping one ONB to another. Finally, we examine 
the matrix for U with respect to the orthogonal decomposition given 
by the Cuntz operators. 

3.1. The Cuntz relations. It is known that for every integer N > 1, 
there is a rich variety of non-equivalent representations of the Cuntz 
relations on N generators |Cun77] . The corresponding C*-algebra on 
the N generators is denoted On- Even the case N = 2 offers a rich 
family of such representations with applications to analysis of wavelets 
and, more generally, to fractals generated by two transformations; see 
Definition 13. II below. We note that the use of the Cuntz relations is not 
restricted to the harmonic analysis of affine fractals; the relations can 
be applied much more generally. Our purpose here is to isolate a par- 
ticular representation of O2 which lets us produce powerful algorithms 
for Fourier expansions in L 2 (n). 

Definition 3.1. We say that operators So, Si on L 2 (fi) satisfy Cuntz 
relations if 

(1) S S* + SiS* = J, 

(2) S*S j = 5 iJ I iorij = 0,1. 
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When these relations hold, we have a representation of the Cuntz al- 
gebra 2 . 

We begin by defining scaling operators So and Si on the Hilbert space 
L 2 (/i) and showing that they satisfy Definition 13. 1( this representation 
yields an orthogonal decomposition of L 2 (fi). To motivate the con- 
nection between the harmonic analysis of afline fractals and operator 
representations within a C*-algebra, we look to the recurring notion 
of self-similarity. The self-similarity or scaling invariance present in T 
given by 

r = 2nTU (2nT + ^y 

has a representation via the Cuntz relations on bounded operators on 
L 2 (/j,). The introduction of representations of the Cuntz relations al- 
lows us to transfer spectral duality questions to a framework in geom- 
etry of Hilbert space. 

Given the Bernoulli measure \i with scale factor let T be the 
canonical spectrum. We define So, Si on the canonical ONB E(T): 

(3.1) S e 7 
and 

(3.2) 5ie 7 = e™ + 2n 7 - 

Using Parseval's identity, we find that So and Si are isometries since 
they each map E(T) into a subset of itself. The range of So is the 
subspace spanned by E(2riT), and the range of Si is the subspace 
spanned by E(^ + 2nT). It is apparent that the range spaces are 
orthogonal and that L 2 (/j,) is the direct sum of these spaces. We can 
extend this decomposition by observing that So and Si satisfy the 
Cuntz relations from Definition 13.11 

In order to verify that So and Si verify the Cuntz relations, we first 
compute their adjoints. Given £, 7 from T, 

(e s , So*e 7 ) = (S e € , e 7 ) = £(2n£ - 7). 

If 7 is not in 2nT then the exponentials are orthogonal. If, on the other 
hand, 7 = 2n7' for some 7' G T, then by Equation (12. 7p . 

£(2n£ - 2717') = //(£ - 7') = (e ? , e y ). 

Thus, we have 



(3.3) S e 7 



csl when 7 G 2nT 



2n 



otherwise. 



s 
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Similarly, 

(e £ , S*e 7 ) = (5ie s , e 7 ) = Jl(2n£ + - - 7) . 
Here, we get zero when 7 e 2nl\ If 7 = 2n7' + |, 

(l(2nZ + -- 2n T ' - -J = ju(2n(£ - 7 ')) = - Y) = fe, ey). 
This gives 

(3.4) S*e 7 



e 7 _« when 7 £ | + 2nr 

2n 

otherwise. 



Proposition 3.2. The operators So and Si defined in Equations A3.1\) 
and $3.2) satisfy the Cuntz relations in Definition \3.1\ 

Proof. We start with the second part: S*Sj = 5ijl for i, j — 0, 1. Then 
we show that So^Sq + SiS* = I. 

(2) It follows directly from the computation of the adjoints Sq and 
SI that applying either Sq to an element in the range of S\ or applying 
SI to an element in the range of So yields 0. Also from the definitions, 
we see that the adjoints act as left inverses to the original operators, 
i.e. S*S = I = S^Sx. 

(1) Given the isometries Si, i — 0,1, SiS* is a projection onto the 
range of Si since 

(s l s;)(s l s:) = s l (s*s 1 )s: = s l s* 

by (2) above. We have already shown that the ranges are orthogonal 
and their direct sum is all L 2 (n), hence SqSq + S\S* = I . □ 

Remark 3.3. Proposition 13.21 can be seen as a special case (N = 2) of 
Theorem 2.3 in |DJllj . which uses the Cuntz relations to formulate 
an algorithmic approach for constructing different Fourier ONBs for a 
fixed IFS measure \i. We note here, however, that the spectral scaling 
by p which we introduce in Section H] to build the operator-fractal U is 
not part of a Cuntz algebra representation. 



We will use S and Si to describe an orthogonal decomposition of 
L 2 ([i). The first step will be to define the following decomposition of 
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the canonical spectrum T. Let 

T = r \ 2nY = - + 2nT 

r x = 2nT \ (2n) 2 r = 2n Q + 2nr) 

T fc = (2n) fc r\(2n) fc+1 r = (2n) fc ^ + 2nr). 
It is readily verified that the sets Tfc are disjoint and form a partition 

O fr\{o}. 

The elements of T are in a natural one-to-one correspondence with 
infinite sequences of bits from {0, ^} having finitely many nonzero el- 
ements: 

m 

(3.5) 7 = ^2 a ^ 2n ) 1 ( fl 0' °i> a 2> • • • > a m, 0, 0, 0, . . .), 

where aj G {0, |}. We can describe as the set of all elements from V 
which correspond to a sequence having exactly k leading zeroes (note 
the index starts at 0): 

Tl 

(0, 0, ... , 0, -, a k+1 , . . . , a m , 0, 0, 0, . . .). 

The operators So and Si can also be associated to operators on 
sequences via their behavior on the ONB elements e 7 . The operator 

50 takes e 7 to e2n 7 , which corresponds to a right shift on the sequence 
associated to 7. Since Si takes e 7 to e™ + 2n 7 , we associate the action of 

51 to the operator that shifts the sequence for 7 right and places | in 
the vacant position. 

7 -H- (a ,ai,a 2 , ■ ■ ■ ,a m , 0,0,0, . . .) 
2^7 ^ (0, a , ai, . . . , a m , 0, 0, . . .) 

Tl Tl 

- + 2n7 o (-,a ,ai,...,a m ,0, 0, ...) 
We combine the associations above to prove the following. 

Lemma 3.4. The closed span of E(Tk) is the subspace SqSi(L 2 (/j,)), 
for each k — 0, 1, 2, . . .. 

Proof. We show that G E(Tk) if and only if eg = S^Sie^ for some 

7 g r. 

From the definition, 7' G Tfc if and only if 7' = (2n) fc (| + 2717) for 
some 7 G T. Then, for this 7, 

5'oS'ie 7 = Soeri + 2rvy = e (2n) fc (f +2n 7 ) = <2 7 '. 
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Thus, 

E(T k ) = 5 *5 1 (£(r)), k = 0,1,2,.... 

Since E(F) is an orthonormal basis for L 2 (fi), the closed span of E(Fk) 
is the subspace 5q5i(L 2 (/i)). □ 

For k — 0, 1, 2, . . ., let W 7 ^ denote the following set: 

W fe :=5*5 1 (L 2 (//)). 

Note that Wk+i = SoWk- By the previous lemma, each Wk is spanned 
by E(Tk). Since the sets form a partition of T\{0}, the Wk subspaces 
are mutually orthogonal, and their direct sum yields L 2 (/i)0sp{e o }. We 
record the orthogonal decomposition of L 2 (fi) in the following propo- 
sition. 

Proposition 3.5. Given the spaces 

(3.6) W k = 5 fe 5!(L 2 (/i)) for k = 0, 1, 2, ... , 
we have 

oo 

(3.7) L 2 (n) = sp{e }®Q)W k . 

4. The operator-fractal U 

There are many cases (for example, [LW021 lDJ09l IDHS091 IJKSllj ) 
in which a fractal L 2 (fj,) Hilbert space has more Fourier bases than just 
the canonical one. We will refer to these as alternate ONBs for the 
space. One technique in the papers above for finding alternate ONBs 
is to scale the spectrum V of the canonical ONB by an odd integer p. 
The new set pT always yields an orthogonal collection of exponentials 
and often yields an alternate ONB. Examples include p = 5 for the 
A = \ measure [DJ091 Prop. 5.1] and p = 3 for the A = | measure 
j.TKSlll Thm. 3.5]. 

In this section, we will use pairs of canonical and alternate ONBs 
to define the unitary operator U. We will demonstrate that U has 
inherent fractal-like properties, and we will find additional structure in 
the A = j, p — 5 example. 

Consider an arbitrary pair for A = where the canonical set Y = 
T(j-) and the scaled set pT are both spectra. Define the operator U 
on T by 

(4.1) Ue 1 — e pl . 

Since U maps an ONB to another ONB, it is a unitary operator on 
L 2 (fi). Note that ?7 2 (e 7 ) ^ e p 2 ry in general. Computation of powers 
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of U requires an infinite expansion with respect to the E(T) ONB, 
i.e. U 2 e 1 = J^teri^'Ue^Ue^. Toward this end, it becomes useful to 
study the matrix representation for U with respect to E(T). We find 
a surprising self-similarity in the matrix for U. 

Given 7, £ G T, the matrix entries for the operator U are 

^(6,7) = (ef.^B,) =^(P7~0- 

If we reorder the index set T by the subsets T^, we find that the matrix 
for U has a block structure. The block structure arises because the 
subspaces Wk and U{Wk>) are orthogonal for k ^ k', but we will in fact 
prove the stronger result that each space Wk is invariant under U. 

4.1. Invariance properties of U. In this section, we examine the 
interplay of U and the Cuntz subspaces. We find that the subspaces 
are invariant under U. From now on, we use £ to denote 

00 

C = L 2 (fi)Qsp{eo} = ^W k . 

Lemma 4.1. For each k > 0, Sq(C) is invariant under U . 
Proof. 

Case 1: k — 0. We show that Ue^ is in the orthogonal complement of 
sp{eo} if and only if 7 7^ 0. This guarantees that Ue 1 G C 

If 7 = then t/e 7 = eo- Conversely, if 7 7^ 0, then Ue 1 — e P7 where 
P7 is a nonzero element of pT. Then, since E(pT) is an ONB which 
contains e , we have (e P7 , e ) =0. 

Case 2: fc > 1. Let 7 G T \ {0}. By the definition of S in dXTJ and 
by Equation f)2.5p . 

= ]T /I(p(2n) fc 7-ne 5 '- 
«'er\{o} 

In the expansion above, we leave out £' = since USqC C USqC C £. 

Consider the terms in the sum for which £' ^ (2n) fc r. Recalling 
that the sets 1^,1 G N partition T \ {0}, we find that £' £ (2n) k T 
is equivalent to £' being an element from r io for some %q < k. Then 
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£' = (2n) io (f + 2n£) for some (6 T, and 

p(2n) fc 7-r = ^(p(2n) fc 4 7 -4(2nr(^ + 2<)) 

= - (2n) io (p(2n) fc " io 4 7 - 2n - 8n^j 

= J (2n) io+1 (p(2n) fe ^°- 1 4 7 - 1 - 4f) 

G 2, the ju zero set from (12.61) . 

Recall that, whether n is even or odd, 4 7 and 4£ are even integers. 
Therefore, the quantity in parentheses in the last expression of the 
equation above is an odd integer. 

This proves that /K>(2n) fc 7-£') = when f £ (2n) k T, hence US k e y 
is in the closed span of E((2n) k (T \ {0}), which is exactly the space 
Sq(C). Therefore Sq(C) is invariant under U . □ 

Next, we show that the Wk subspaces are invariant under U. 

Proposition 4.2. Each space Wk = S'o>S'i(L 2 (/x)) is invariant under 
U. 

Proof. Let ey G Wk- This implies that ey G Sq(jC) C £. We showed 
in Lemma E31 that = sp{E(Tk)}, so we have 7 ' G is of the form 
7 ' = (2n) fc (^ + 2n 7 ) for some 7 G T. Then, 

f/ey = e p y = ^ £(pY - £')e € /, 
£'er\{o} 

where, by Lemma [4. II the coefficients in this sum are possibly nonzero 
only for nonzero £' G (2n) fc r. We will show that the coefficients are 
also zero for £' G (2n) fc+1 T C (2n) fc r, hence the only possibly nonzero 
coefficients are from £' G IV 

Let f G (2n) fc+1 r, i.e. f = (2n) k+1 £ for some ( G T. Then, e e G 
S +1 (£) and 

W -O = £(p(2")*(£ + 2n 7 ) - (2n) fc+1 e). 
We then compute that 

p(2n) fe (- + 2n 7 J - (2n) k+1 £ = -(2n) k+l (p + 4p 7 - 4f) G 2. 

Since pj' — £' is in the zero set, Uey is orthogonal to each e^> where 
£' G (2n) fc+1 r. Thus, we have C/ey G IV^, so Wk is invariant under 
[/. □ 
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4.2. Commutation relations for U and Sq. We can actually say 
something stronger than invariance about the interaction between U 
and So- When n is even, the operators commute. When n is odd, they 
don't quite commute but do have a straightforward relationship. We 
note that U and Si do not commute. 

Theorem 4.3. If n is even, then the operators Sq and U commute on 
L 2 (/i). 

Proof. Let 7 G V . Note first that 2n7 G T, while py is not necessarily 
an element of T. We must therefore expand both expressions USqc^ 
and SqUc 1 in terms of the ONB E(Y). First, 

(4.2) " 6r " er 
r?er 

On the other hand, 

(4.3) USoe^ = Ue 2ni = e p2nT 

Since T = 2nT U ^| + 2rarj = 2nT U To, we can split the sum over T 
into two terms: 

U S e^ = ^ KP 2n l ~ e ? 

(4.4) ?er 

= Yl KP 2n l - e 5 + KP 2n l - e ?- 
£e2nr £er 



However, by Lemma [4.11 Ji(p2n'y — £) = when £ G r , hence 

U S e 7 = Y HP 2n 7 ~ 

(4.5) ?€2nr 

e'er 

Since n is even, 7 G Z and £ G Z. By applying Equation ( 12. 7p . we have 

(4.6) U S e y = Y A*(P7 ~ O e 2<' , 

which agrees exactly with the expression for S^Ue^ in Equation ( 14. 2p . 
Note that U and So commute when restricted to £ as well. □ 

In the previous theorem, we depended on the fact that n was even 
to apply Equation (12 .7p in the last step in Equation (14.51) . In fact, if 
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n is not even, U and Sq do not quite commute, as the next theorem 
shows. First, we point out some straightforward observations. 

• When n is odd and 7 6 T = | + 2nT, we see that 7 is not an 
integer, and is in fact an element of |(2N — 1). Therefore, 2rry 
is an odd integer for 7 G Tq and an even integer for 7 G T \ Tq = 
2nT. 

• Let n and p be odd and let 7 G To and ( G T. Then, since 
2t2 7 e N, 

P7-2< G -(2Z+1). 



Given n and p odd, 



cos(2tt(p7 - {)) 



1 7, f G 2nr 

-1 7 e2nr, eer 
[-1 7 er , ee2nr. 



Theorem 4.4. Suppose n is odd and 7 G T. Write 

SoUe^ = s + w, 
where s G Sq(L 2 (/i)) and w G W\ — S S 1 (L 2 (fi)). Then 



(4.7) 



U S n e. 



e 7 



s — w i/7 6 2nT 

-s + w z/ 7 g f + 2nr = r . 



Proof. By Equation ( 14. 2 p in Theorem 14. 3[ we have 
(4.8) So?/ e 7 = ^2 fiivy ~ v) e 2n V 

Since 2nT = (2n) 2 T U r 1; we split up the sum into these two pieces, 
which we call respectively s and w. 



SqU e 7 = ^2 fi(p ■ 1 — 2n ■ rj 



??er 



(4.9) 



+ Yl ' 7 ~ \ ~ 2n ' y e 2n(f +2nr,) 

r?er 

s®w G 5 2 (L 2 (/j)) © Wi. 
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We split the sum for USoe y in the same fashion, yielding s and w: 

U 5 e 7 = ^ V>(P 2n l - e ? 

£€2nr 

(4.10) ?er 



^2 / r 2 , 



s©w G S$(L 2 (ji)) ®W : 



Case 1: 7 G 2nr. Let 7 = 2n7' G 2nT. Then, making use of Equation 
(12.71) . we see that the coefficients for s in Equation (I4.10p are 

/2(p(2n) V - (2n) 2 
= cos (2^(^2727' — 2n£)) ji{p2n^' — 2n£) 
= Ji{p2n^' — 2n£), 

since p2n r )' — 2n£ G Z. These are the same as the coefficients for s 
from Equation (14.91) . so s = s. 

The coefficients for w in Equation (14.101) are 



ji(p(2n) 2 i -2n(^ + 2n^ 
= cos(2n(p2n-f' - | - 2nf)) Ji(p2rii - | - 2nf) 
= (-L)/2(p2n T ' - ^ - 2n^ 



The cosine factors are all —1 since p2n^' — | — 2n£ is an element of 
|(2Z + 1). Therefore, the coefficients for w are exactly the opposite 
sign from those for w, i.e. w = —w. 
Case 2: 7 G T = § + 2nr 

Let 7 = f + 2717'. Again using Equation (12. 7ft . we find that the 
coefficients for Is are 



$(p2nQ + 2n^ - (2n) 2 ^ 
= cos^27r(| + 2ny - 2<jj + 2n l' ~ 2n £ 
= (-l)/2(^ + 2n 7 '-2<) 
The coefficients for s are the opposite sign to those for s, so 7s = —s. 
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The coefficients for w are 



#(p2nQ + 2nY) - 2n(| + 2< 

= (2n ^^-^ — n + p2n^' — 2n£ 

^/p — 1 
= /i^ — - — n + p2ny — 2n£ 



since is an integer. These are the same coefficients as those for w, 
so w = w. □ 

Corollary 4.5. When n is odd, the commutator USq — SqII is given 
by 

-2w 7 G 2nr 



(^o - S U)e, 



-2s 7 e r c 



where s and w are as defined in Proposition \4-4 



The isometry operators Sq, Si have range projections SqSq and S\SZ 
respectively. The preceding results lead us to find that U commutes 
with the range projections. 

Corollary 4.6. The unitary operator U commutes with SqSq and S\S*. 

Proof. This follows immediately because the spaces Wo and So(L 2 ([i)) 
are invariant under U. □ 

4.3. The matrix for U. We showed in Proposition 13.51 that the Wk 
subspaces are orthogonal, and since Proposition 14.21 showed that each 
Wk is preserved by the operator U, we find that the infinite matrix for 
U, when T is ordered by the Tk sets, has a block matrix structure. We 
note that, with the exception of the upper left corner, all of the blocks 
are infinite. 








r 




r 2 


r 3 







i 
















To 





U Wo 













Ti 








u\ Wl 










r 2 











U\w 2 







r 3 














U\ Ws 
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Next, we show that, in fact, the infinite blocks above are all the 
same. For each k G N , consider the matrix of the operator U\w k with 
respect to the indexing set T k . Let £',7' G T k , so they have the form 



£' = (2n)*Q + 2f* 
y = (2n) fc Q + 2n 7 



and 



where £,7 G T. In the next lemma, we show that each (£',7') entry 
of U\w k is equal to the corresponding (£',7') entry for U\w - In other 
words, the infinite blocks in the matrix for U are all the same. 

Lemma 4.7. For all k G No, the infinite matrix for U\w k has the same 
entries as the matrix for U\w - 

Proof. Define £',7' as above. The (£',7') entry of the matrix for U\w k 
is 

U\w k {£,l') = (e(2n) fc (n/2+20'^ e (2n) fe (n/2+2n 7 )) „, ■ 
\ / L 2 (/i) 

We now compute U\w k (£',j')' 

e (2n) t (n/2+2n()) U & (2n) k {n / 2+2n~f) ) „ 

(4.12) =V>[p- (2n) fc (| + 2n 7 ) - (2n) fc (| + 2< 

= $ (^2n) k ^P~^ n + 2n P1 - 2< 

By the invariance property of \x in (12. 7p . 

ju(2nt) = cos(27rt)ju(t). 

If t is an integer, then fi,(2nt) = fi(t). Since p — 1 is even, the number 
t given by 

(p — l)n 

t = v ' + 2np7 - 2ni 

in Equation ( 14.12ft is a whole number, and thus we can disregard the 
k powers of 2n. Therefore, we reach our desired conclusion: 

U\w k (C',l') = ri[{2n) k ( [ ^^ + 2np 1 -2ni 
fi[ ( p ~^ n + 2np-f - 2< 



u\w (e,i') 
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since, in Wq, £' 



2n£ and 7' = ^ + 2wy. 



□ 



Since the spaces Wk are mutually orthogonal and are each invariant 
under the action of U, the infinite matrix representation for U has the 
following block form. As before, the index sets have been written across 
the top row and left column. 

Theorem 4.8. The matrix representation of the operator U , for the 
given ordering of the index set T, has a block diagonal structure of the 
form 



(4.13) 



U 








r 


Ti 


r 2 


r 3 







1 
















To 





U\w 













Ti 








U w 










r 2 











U Wo 







r 3 














U Wo 



















In other words. 



U ~ P. 



fc=0 



where P eo is the orthogonal projection onto sp{eo}. 



Remark 4.9. Making use again of Equation (12. 7p . and depending on 
the values of n and p, we can write the entries of the matrix blocks as 



I — ^" + 2np-f - 2n£ 



(2n)(^+ P1 - 



p — 1 



4.4. The individual blocks in the p = 5, n = 2 case. In this section, 
we consider the matrix of U in the special case p = 5 and n = 2. Here, 
U{e 1 ) = e 57 , and we zoom in on the subspace W = Si(L 2 (/ii )). We 
will study the structure of one of the blocks in the matrix form for U 
in Equation (I4.13p . We find that, with a specified ordering of the index 
set of a block defined by the Cuntz operators So and 6*1, there is again 
a self-similar block structure in the matrix of {7|w - 
When n = 2, the zero set for ju is 



(2£ + r 



GZiGN, 
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The canonical ONB is 

r=|^a i 4 i : a* e {0,1} J = {0,1, 4, 5, 16, 17, 20, 21,...}. 

By (I4.12p . the (£',7') entries in U\w k are given by 

(4-14) tfM£',Y)=£i(i + 57-0, 

4 

where £' = 4 fc (l + 4£) and 7' = 4 fc (l + 47). Each such block is indexed 
by the associated T k , and by Theorem 14.81 the entries do not depend 
on k. It is sufficient, then, to examine the first block, indexed by IV 

We recall that the entries for the matrix of U are £/(£, 7) = j2(5j— £). 
This yields the expression of the restriction U\w as unitarily equivalent 
to the product of U and a multiplication operator, which brings us to 
our surprising operator-fractal property of U. 

Theorem 4.10. Let M ei be the multiplication operator by the function 
ei on L 2 (fi), and let P eo be the orthogonal projection onto sp{e }. Then 
U has the form 

00 

(4.15) U ~P eo ®Q)M ei U. 

k=0 

Proof. First, observe that 

U\w — S^U Si, 

since the range of the isometry Si is exactly the subspace W . We see 
that 

(M ei C/e 7 ,e e ) = (e 57+ i,e € ) 

= /2(l + 5 7 -0 
= (SlUSie.y,e^) 

- U\ Wo (-y,$,) by Equation (gH]). 
The result is the following commutative diagram. 



£ 2 (/ii) W 

4 



M ei U 



4 

The result then follows from Theorem 14.81 □ 
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The exponentials E(Tq) form an ONB for the subspace Wo = Si(L 2 (/ii )). 
We make use of the Cuntz relations on the operators S and Si to divide 
W into a direct sum of subspaces. Let 

(4.16) V k = S 1 S k S 1 (L 2 (^)), k = 0, 1, 2, . . . . 

4 

Each V k has ONB E(Tf.), where we define 

f fc = 1 + 4 fc+1 (l +4r). 

Using the correspondence from Equation (13. 5ft between elements of T 
and infinite sequences, an element 7 from Tk corresponds to a sequence 
where ao = a k +i = 1 and a± = ■ ■ ■ = a k = 0: 

7 -H- (1, 0, 0, . . . , 0, 1, a k+2 , a k+3 , . . . , a m , 0, 0, . . .). 

We see that T k C T for all fc 6 No and that Tj, T k are pairwise 
mutually disjoint. Moreover, the T k sets form a partition of T \ {1}: 

00 

(4.17) r = {i}uyr fe . 

k=0 

Making use of Equation ( 14.161) . we have the following decomposition of 
W : 

00 

(4.18) W = S^ifi)) = sp{ ei } © V k . 

k=0 

We now compute the entries of the matrix for {/|wv We first observe 
that the majority of the blocks of the matrix are zero, and then we focus 
on the four row and column blocks with possibly nonzero entries. 

Lemma 4.11. Let i,j > 1. Let f e f % = 1 + 4 i+1 (l + 4r) and 
y e fj = 1 + 4^ +1 (l + 4r). Then 

so the (Ti,Tj) block in the matrix of U\w is the zero matrix. 

Proof. Let £' = 1 + 4 m (l + 4f) and 7' = 1 + 4^' +1 (l + 4 7 ) for ^, 7 e Y. 

Since U\w {£,' , 7') = jwi(5y — ^'), we compute 57' — 

4 

5 7 ' - = 5(l + 4 4+1 (l + 4 7 )) - (l + 4 J+1 (1 + 40) 
= 5 + 5 • 4 m (l + 4 7 ) - 1 - 4 j+1 (l + 40 
= 4(l + 4 4 (l + 4 7 ) -4 J '(l + 40) 
= 4 ■ odd integer e Z{jLi). 
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□ 

Lemma 4.12. Let 7' G r = 1 + 4I\ The following hold for the matrix 
entries in the first row ofU\w - 

(a) Block (1,1): U\ Wo (l, 1) = 0; 

(b) Block (l,r ): Ifi GT = 1 + 4(1 + 4r) ; then U\ Wo (l, 7') ± 0; 

(c) Block (1, f fc ), fc > 0: J/Y G V k = 1 + 4 fc+1 (l + 4r) for k > 0, then 
U\ Wo (l,j') = 0. 

Proof. We work with the equation [7(1,7') = filial' — 1) for 7' G IV 

4 

(a) We have 5(1) - 1 = 4 G by (12^]) . so U\ Wo (l, 1) = 0. 

(b) Since 7' G To, we have 7' = 1 + 4(1 + 47) for some 7 G T. Then 

5 7 '- 1 = 5(1 + 4(1 + 4 7 )) -1 

= 4 + 20 + 8O7 = 4(6 + 2O7) 4 Z(a L ). 

4 

(c) We repeat the above computation, but now k > 1: 

5 7 ' - 1 = 5(1 + 4 fc+1 (l + 4 7 )) - 1 
= 4 + 5 • 4 fc+1 + 5 ■ 4 fc+2 7 
= 4(1 + 5 • A k + 5 • 4 fc+1 7 ) 
= 4 • odd integer G Z(pi). 

This proves /2i(57' — 1) = 0. 

4 

□ 

We now turn to the first column of the matrix of £/|w - 

Lemma 4.13. The following hold for the matrix entries in the first 
column ofU\w - 

(a) Block (ITo, 1): This block contains nonzero entries. 

(b) Block (Tk, 1), k > 0: Each of the blocks (I\,, 1) is zero. 

Proof. Let f = 1 + 4 fc+1 (l + 4f ) , where f G T. 
(a) Suppose = 0. We compute 5 — 

5 - (1 + 4(1 + 40) = 4 - 4(1 + 40 
= 4(1 -(1+40) 
= -4 2 e 

If ^ G T is even, which it can be (£ = 4, 16, 20, . . .), the associated 
matrix entry is nonzero. 
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(b) We compute 5 - f = 5 - (1 + 4 m (l + 40), where k > 1: 

5 - (1 + 4 fc+1 (l + 40) = 4 - 4 fc+1 (l + 40 

= 4(l-4 fe (l+40) 

= 4 ■ odd integer £ Z(pi). 

4 

□ 

Finally, we examine the (T^To) and (To,!^) blocks in the matrix 
for Z7]w - 

Lemma 4.14. Lei k > 1. 

(a) Block (r ,r ): 77ms fr/oc£; is zero. 

(b) Block (ITo, r&): 77ms 6Zoc& contains nonzero entries. 

(c) Block (r^jTo): 77ms WocA; contains nonzero entries. 

Proof. Again, we use the formula U\w (£',j') = — £')• 

4 

(a) Let = 1 + 4(1 + 40 and 7' = 1 + 4(1 + 4 7 ) , where f , 7 £ T. Then 

5 7 ' - = 5(l + 4(1 + 4 7 )) - 1 - 4(1 + 40 
= 5 + 5 • 4(1 + 47) - 1 - 4(1 + 40 
= 4(l + 5(l + 47) -1-40 ) 
= 4(5 + 207-4^) eZQjlx). 

(b) Let f = 1 + 4(1 + 40 and 7' = 1 + 4 fc+1 (l + 4 7 ), where f , 7 £ T 
and k > 1. Then 

5 7 ' - = 5 (l + 4 fe+1 (l + 4 7 )) - 1 - 4(1 + 40 
= 5 + 5 • 4 fc+1 (l + 4 7 ) - 1 - 4(1 + 40 
= 4(l + 4 fc (l + 4 7 )-l-4£) 

= 4(4 fc (l + 4 7 )-4£) 

= 4 2 (4*- 1 (l + 4 7 )-e) 

which is not necessarily a zero of jui. For example, if fc and 7 

are fixed, then 4 fc " 1 (l + 47) £ T. Set £ = 4 fc-1 (l + 47). Then 
57' — £' = (), which does not belong to the zero set of jt/i. 
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(c) Let f = 1 + 4 fe+1 (l + 40 and V = 1 + 4(1 + 4 7 ), where (,7 G T 
and fc > 1. Then 

5 7 ' - £' = 5 (l + 4(1 + 4 7 )) - 1 - 4 fc+1 (l + 40 
= 5 + 5 ■ 4(1 + 4 7 ) - 1 - 4 fc+1 (l + 40 
= 4(l + 5(l+4 7 )-4 fc (l + 40) 

= 4(6 + 20 7 -4 fc (l + 4O), 

which cannot belong to the zero set of jE/i; 2(57' — 0) is congruent 
to 2 (mod 4). 

□ 

The results from our previous lemmas combine to give a block struc- 
ture to the matrix for the operator U\w - 

Theorem 4.15. The matrix for U\w has the following block form, 
where the index set To for the ONB of Wo is ordered by {1} U [JlV 
The notation * indicates an entry or block that is not necessarily zero 
and indicates an entry or block that contains all zeros. 



(4.19) U\ Wo ~ 





1 


r 




r 2 


r 3 




1 





* 













r 


* 





* 


* 


* 









* 













r 2 





* 













r 3 





* 
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